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Abstract 



The Debye mass mn is computed nonperturbatively in the deconfined phase 
of QCD, where chromomagnetic confinement is known to be present. The latter 

cn \ defines m£) to be niD = cd^/o~s^ where cd — 2.06 and as = crs{T) is the spatial 

string tension. The resulting magnitude of rao[T) and temperature dependence 

f^ , are in good agreement with lattice calculations. Background perturbation theory 

^D \ expansion for moiT) is discussed in comparison to standard perturbative results 

^I ■ and recent gauge- invariant definitions. 

o 
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"Fi : 1 Introduction 

>. 

^ ■ The screening of electric fields in QCD was originally considered in analogy to QED 

j^ . plasma, where the Debye screening mass was well understood PP, and the perturbative 

leading order (LO) result for QCD was obtained long ago [2 , m\) = ( "^ + -g^ ) qT- 
For not very large T, however, the purely perturbative expansion is not reliable, and 
attempts have been made to use the effective "id theory jHl IH El El to define the Debye 
mass m£, through the coefficients, which are to be determined nonperturbatively |2|- In 
doing so one obtains a series [TJ, with the leading term of the same form as m]-, . 

The lattice calculations of rriEiiT) have been made repeatedly |H]-|IS], and recently 
mniT) was computed on the lattice for A^/ = 0, 2 [T31II1| using the free-energy asymptotics 

<5Fi(r,T) = Fi(r,T) -Fi(oo,T) ^ _'^'^^^-mn(T)r ^ (^^ 
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where Fi{r,T) was found from color singlet Polyakov loop correlator. A comparison of 
lattice defined mniT) with m^ made in [15^ in the interval from Tc up to temperatures 
about 5.5Tc shows that one requires a multiplicative coefficient Ajq^^Q = 1.51, Aiy^=2 = 
1.42. A difficulty of the perturbative approach is that the gauge-invariant definition of 
the one-gluon Debye mass is not available. The purpose of our paper is to provide a 
gauge-invariant and a nonperturbative method, which allows to obtain Debye masses 
in a rather simple analytic calculational scheme. In what follows we use the basically 
nonperturbative approach of Field Correlator Method (FCM) |TBj-|221 ^^^ Background 
Perturbation Theory (BPTh) for nonzero T [231211123 to calculate m£,{T) in a series, 
where the first and dominant term is purely nonperturbative, 

mD{T) = Mo + BPTh series. (2) 

Here Mq is the gluelump mass due to chromomagnetic confinement in 3d, which is 
computed to be Mq = cd^/o~s, with as{T) being the spatial string tension and cd — 2.06 
for A^c = 3. The latter is simply expressed in FCM through chromomagnetic correlator 
|19j . and can be found either from lattice measurements of the correlator itself as in [TTj . 
or from the 3d effective theory |3]-|Zj, \f^s = c„g'^{T)T, or else from the lattice data pBj . 
Therefore Mq(T) is predicted for all T and can be compared with lattice data jISl Cll 113 , 
see Fig. |21 

We note, that m£)(T) is defined here as the screening mass in the static QQ potential 
Vi, which can be expressed through the gauge- invariant correlator of chromomagnetic and 
chromoelectric fields [211 123 I2H] • The screened Coulomb part of the potential Vi coincides 
with the singlet free energy Fi{r, T) at the leading order [22], and in what follows we shall 
consider also the leading order in BPTh, where the static potential Vi{r) has a term of 
the same form as the r.h.s. of Eq. (^. 

The paper is organized as follows. In section 2 the nonperturbative part and the per- 
turbative BPTh series for the thermal Wilson loop are defined, and the gluelump Greens 
function is identified, using the path-integral formalism. In section 3 an effective Hamil- 
tonian is derived and the first terms of expansion ^ are obtained for m£,{T) computed 
through the spatial string tension o"s(T). In section 4 a comparison is made of mij{T) 
with lattice data and other approaches. Section 5 is devoted to a short summary of results 
and outlook. 



2 Background Perturbation Theory for the thermal 
Wilson loop 

It is well known that the introduction of the temperature for the quantum field system in 
thermodynamic equilibrium is equivalent to compactification along the euclidean "time" 
component X4 with the radius j3 = 1/T and imposing the periodic boundary conditions 
(PBC) for boson fields (anti-periodic for fermion ones). Thermal vacuum averages are 
defined in a standard way 

{...) = ^f [DA]...e-'^^^\ (3) 

^/3 JPBC 



where partition function is 

Z^= /" [DA]e~^f'^^\ Sp= f dxi I (fxLYM. (4) 

JPBC Jo J 

One starts as in |28j and |30| with the correlator of Polyakov loops (L(0)L+(r)) 
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L(x) = -^trPexp(z(7 / y44(x, X4)(ix4), 
and obtains (cf. 



/3 
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(L(0)L+(r)) = — exp (-/?Fi(r, /5)) + ^^^ exp (-/3F8(r, /?)) . (5) 

As it is explained in the Appendix in [2Z| the representation (0) can be obtained 
from two Polyakov loops by identical deformation of contours with tentackles meeting at 
some intermediate point and subsequent merging of contour into one Wilson loop using 
completeness relation at the meeting point 5a^i3^5a2i32 = -^^01/32^02^1 '^'^'^%oJ''p-i_02^ ^^^ ^^^^ 
term contributing to the free energy Fi of the static QQ-^gaXi in the singlet color state, 
the second to the octet free energy Fg. Accordingly one ends for Fi with the thermal 
Wilson loop of time extension (3 = 1/T and space extension r, 

exp(-/3Fi(r,/5)) = {W{r,P)) = ^JtiPexpitg J A,dz,)). (6) 

Note, that in contrast to the case of the zero-temperature Wilson loop, the averaging in 
© is done with PBC applied to A^, as in ©, Q. 

Eq.® is the basis of our approach. In what follows we shall calculate however not Fi, 
which contains all tower of excited states over the ground state of heavy quarks QQ, but 
rather the static potential Vi(r, T), corresponding to this ground state, for more details 
see EH- 

Separating, as in BPTh |23^ the field A^ into NP background 5^ and valence gluon 
field a^, 

A^ = B^ + a^ (7) 

one can assign gauge transformations as follows 

B^^U+{B^ + -d^)U, a^^U+a^U. (8) 

g 

As a next step one inserts ((Tj) into (jH)) and expands in powers of ga^, which gives 

{W{r,P)) = {W^'\r,P))B + {W^'\r,P))B,a + • • • , (9) 

where according to J2SI one can write {T)a = ((r)a)B, and {W^'^^) can be written as 

{W^'^)B,a = ^ [ {tTP<l>([[){a,{x)a,{y))Ml[))Bdx,dy.. (10) 



Here $(11) ^^^ "^(U) ^^^ parallel transporters along the pieces of the original Wilson 
loop W{r,i3), which result from the dissection of the Wilson loop at points x and y, see 
FigHJ Thus the Wilson loop W^'^\r,P) is the standard loop W^^\r,P) augmented by the 
adjoint line connecting points x and y. It is easy to see using (jHl), that this construction 
is gauge invariant. 




Figure 1: The gluon trajectory (wavy line) and the adjoint surface S"^ (dark region) 
attached to the thermal Wilson loop. 

For OGE propagator one can write the path integral Fock-Feynman-Schwinger (FFS) 
representation for nonzero T as in [221 

Gf,uix,y) = {af,{x)a^{y))a = 
= r ds f {D'zr^yeM-KW'KC.y) (PpeM'^ig f F,p{z{r))dr)] , (11) 
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where the open contour C^y runs along the integration path in pi|) from the point x to 

■^Ixy 



the point y as shown Fig^ and K = ^ J^{z^)'^dT. The path integration measure {D^ 
is given by 



fe=l ^ J J \ J n=-oo \ \k=l 



(12) 
with Ne = s and Azfj^{k) = z^{k) — z^{k — 1). Thus, {D^z)^y is a path integration with 
boundary conditions z^(r = 0) = a;^ and z^{t = s) = y^j, (this is marked by the subscript 
xy) and with all possible windings in the Euclidean temporal direction (this is marked by 
the superscript w). 

We must now average over B^ the geometrical construction obtained by inserting pi|) 
into (Uni), i.e. 

xy 

{^([{W^{C.^ym\]))B = {W^y{r,P))B. (13) 

xy 



One can apply to (jlUj) the nonabelian Stokes theorem, and to this end one has to fix 
the surface bounded by the rectangular r, (3 with the adjoint line passing on the surface. 
The standard prescription of the minimal surface valid for the fixed boundary contours, 
in our case when chromoelectric confinement is missing and only spatial projections of 
the surface enter, leads to the deformation of the original plane surface due to gluon 
propagation, consisting of this original surface plus the additional adjoint surface S^i 
connecting gluon trajectory with its projection on the plane {r,P), see FigU where this 
projection is simplified to be the straight line. The nonabelian Stokes theorem yields the 
area law jinillll for distances r S> A^, A^— gluon correlation length, A^, ~ 0.2 fm 

{W,y{r, (3))b = exp(-a^ Vne) exp{-a^,^Sf, - a^,^Sf,), (14) 

where S ,' are projections of gluon-deformed piece of surface 5*^; into time-like, space-like 
surfaces respectively. 

For T > Tc one has a^ = and one obtains exactly the form containing the gluelump 
Green's function 

{W^'^)B,a = {igYC2if) f dx, [ dy,Gu{r,h), (15) 

Jo Jo 

where ^4 = X4 — y^, C2{f) = {N^ — l)/2Nc and G44(r, t4) is 

G44(r, U) = j^ dsj (D'z)-^ exp(-ir) exp(-af,^.5j). (16) 

In (jl6|) we have neglected the last exponent on the r.h.s. of (jll|) . which produces 
spin- dependent terms found small in |3Tj, for more discussion see Appendix 4 of ^ 33 . 

Thus the gluon Green's function in the confined phase becomes a gluelump Green's 
function, where the adjoint source trajectory is the projection of the gluon trajectory on 
the Wilson loop plane. 

Now in the deconfined phase, T > Tc, where, magnetic confinement takes place in 
spatial coordinates, so that one can factorize as follows {G^u{x,y) = 6^uG{x,y)) 



j^ ds JiDz,)UDzsUy,G^'\0, 0; s) exp (^-^ ^ 



Gix, y) = J^ ds j iDz,)UDzsUy,G^'KO, 0; s) exp ^-- J^ (zi + zi)dT j , (17) 

where G^'^^ (0, 0; s) is the 2d Green's function with s playing the role of time and interaction 
given by the area law term, exp(— a^^S'^). Here 5*^ is the Nambu-Goto expression 



qH 



dr / d/3^wfw'^ - {wiw'i)'^ 
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and w, = z,(r)/3, z = l,2, w', = ^, Wi = ^. 

In (fTTj) one can specify coordinates in such a way, that X4 = 0, 7/4 = ^4, X3 = 0, 1/3 = r 
and xi,2 = yi,2- 



For G'(^) one can write 

G(2)(0, 0; s) = J{Dz,U{Dz,)oo exp (^-i J\zf + zDdr - a^,^sA . (19) 

The path integral (fT^ can be expressed through the Hamiltonian H^'^\ which is obtained 
from the Euchdean action 

A = j^^ dTC{z,, i,) = \ j\zl + zl)dT + a^,^Sf, (20) 

G^'^\x,y-s) = {x\exp{-H^'^h)\y) (21) 

It is easy to derive, that G^'^\0, 0; s) = G^'^\s) behaves at small and large s as 

G(2)(s^0)oc-^; G^^\s^oo)ocM^exp{-M^s), (22) 

where Mq is the lowest mass eigenvalue of H^'^^ . As an explicit example one can consider 
a rather realistic case when interaction term c^^jS^i in (fTU|) is replaced by the oscillator 
term 

and one obtains 

One can see that asymptotic behaviour (|^ is satisfied provided that u = Mq. On 
the other hand, the eigenvalues of iJ*-^-* (when the role of time is played by s), M^ can 
be expressed through M^, where M„ are eigenvalues of gluelump Hamiltonian H^'^^ when 
Eucledian time evolution is chosen along z^. Those will be found in the next section. Since 
dr = dzs/2fx, dz^ = dt and M„ = 2/i„, the Hamiltonian J H^'^^dr = J(H^'^y2n)dz3 = 
J H^'^^dz^ and one has the equality M^ = M^ with the accuracy of ~ 5% known for the 
einbein technic calculations fI^ . 

Calculating {Dz3)or one has 



A similar calculation with {Dzi)]!^^^ yields 
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(24) 



{Dz,)X exp ( -- / zidr ) = ^= > , exp ( -^^^7^ ) , (25) 



and combining all terms one has 



iid^ - 
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1 V^ „,.„/ (t4 + n/5)2' 
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where we have defined G'^^\s) = 4:nsG''^\0,0; s) so that G(2)(s ^ 0) ^ 1, G^^\s -^ 
oo) -^ exp(— Mqs). 

We are now in the position to obtain the screened static (color Coulomb) potential. 
Indeed, identifying in the lowest order in Oi^g"^) in (0), (US)), (fTBj) . (PU]) . one has 
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F,{r,P) = V^'^ir^P) = -CUW / . dt,G{r,U), (27) 

which can be rewritten using (J26|) as 

vi'\r,P) = -^I^/°°$-P ( -'- ] xis,P)G^%s), (2J 

where xl-s, /?) 



X(., /5) = f] dU £ exp (-^^^^) . (29) 

2 n=— OO ^ ' 

Now for large (3 (small T), f3 ^ r, j3Mq ^ 1, one can keep in the sum (|29|) only the term 
ra = 0, which yields Xn=o{s, P) = V4vrs. From (j^ one then can conclude that s ~ r^, and 
for r'^M^ ~ sMq > 1 one can replace G^'^\s) by the asymptotics, &'^^{s) = exp(— Mgs), 
which yields 

vl'\r,T) = -^^^^e-^'°', rT<l. (30) 

r 

In the opposite limit of small /3 (large T), j3 '^ r, one can use the following relation jHS] 
for the sum in (j^I?|) 

n=-oo ^ /A- fc=_oo \ ^' A- / 

which yields for x(s,/?), 

X(s, /?) = 74^^ (l + 0(e~^"'^'^/^')) (32) 

and hence the screened color Coulomb potential V/ (r, T) has the form pn|) also at 
rT 3> 1. We shall assume accordingly that (jHUj) holds for all temperatures and distances 
r > Ag in the order 0{as), and the next section will be devoted to the calculation of Mq. 

3 Nonperturbative Debye mass 

As it was argued in the previous section, the screened gluon propagator is actually the 
gluelump Green's function, defined in (fTBj) . In this section we shall calculate the gluelump 
spectrum and hence the set of Debye masses. This problem is similar to the calculation 
of the so-called meson and glueball screening masses, which was done analytically in |22j , 
and in our present case we must compute the gluelump screening masses. Below we shall 



heavily use the glueball calculation of 25^, simplifying it to the case, when one of the 
gluon masses is going to infinity-thus yielding a gluelump. 

We note, that the role of time is played by the coordinate ^3, (when the third axis 
passes through the positions of Q and Q) . 

So we write z^ = ^3, < ^3 < r, and define transverse vector z^ = (zi, Z2) and ,24(^3). 
Introducing the einbein variable /i [24,, one has 



J 1 /"r 

^ = 2/i, 0<r<s; K = -j dhii{t^){l + ±1 + zl) (33) 

and G(x, y) acquires the form 

G{x,y) = j Dfi jiDh^MDz,):yexpi-A), (34) 



where the action is 

A — L^ \ , 

adj '-'gl 

Proceeding as in (22| one arrives to the effective Hamiltonian representation 



^ = ^ + ^f*55 (35) 



G'(x,?/) = (x|^exp(-i:f„r)|?/), (36) 

n 

with the temperature-dependent Hamiltonian (r^ = |z_l|) 



Hn = VpI + (2^^^)' + ^ad,r±. (37) 

The spatial gluelump masses are to be found from the eigenvalues of the equation 

Hnvt^ = mP^^I (38) 

and for n = the Hamiltonian ()37p has the form 

^o = Vpi + <r^ (39) 

or in the form with einbein variables which will be useful for discussion 

.2 ,, „2 ,, n-H 2»,2 



^ + 2 + ^adjr± - ^ + 2 + ^^ + 2- 



^cmb ^ ^ + ^ + ^H = ^ + ^ + ^^_^ + _ (40) 



The OGE potential, AV = —Saf^/r, will be considered as the small correction. Note 
the difference between two-dimensional distance r±_ entering in the spatial protection of 
the area in the gluelump Wilson loop, S*^, and the 3d distance r entering in the 3d color 
Coulomb interaction in AV. The eigenvalue of ()4()j) . Hq^^^^p = eoip, with af^ = is 

£o{fi,i^) = ^^ + ^= 41 

2 JJw 



and the minimization in /x, v implied in the einbein formahsm j33] yields 

£o(/xo, z/o) = 2^/^ = 3v/^, (42) 

where Os is the fundamental spatial string tension and a^^- = (9/4)o"s for SUi^^. 

One can compare this value with more exact one, obtained from solution of the differ- 
ential equation in pn|) and to this end one can use the eigenvalue of the screening glueball 
mass found in j23] , (which is larger by a factor of -\/2 than that of our gluelump mass, cf . 
Eq. pij) of J2S] and our Eq. (jHH))- In this way one obtains 

£o = 2.82^5^ (43) 

which differs from ()42|) by 6%. 

In the next approximation the OGE potential for the gluelump comes into play. Here 
one should take into account that the gluon-gluon OGE interaction acquires a large NLO 
correction, which strongly reduces the LO result as it is seen in the BFKL calculation 
(see discussion in j2Sl)! and therefore the effective value of a^^ is smaller than in the 
QQ interaction. Specifically, in the gluelump mass calculation at T = [25 the mass of 
the lowest gluelump for d^^ = is M = 1.4 GeV, and it decreases to M Ri 1 GeV, when 
a^^ = 0.15. This latter value of M is in agreement with lattice correlator calculations pT] : 
the same situation takes place in the glueball mass calculation JJ5|, where also a'^^ = 0.15 
and we shall adopt it in our Eq. (jlUj) . The correction of Eq due to a'f in the lowest order is 
easily computed using PUJ): as a result one has A^o = ~(9/^A")«s*^^/oT ~ —B.OSa'fy/d^. 
As a final result we write the Debye mass (lowest gluelump mass Mq = mu) for af^ = 0.15 

mD = eo + Aeo = (2.82 - 5.08af )v^ = 2.06v^ (44) 

4 Numerical results and discussion 

One can now compare our prediction for mj:,{T) = co ^Jas (T) with the latest lattice data 
J15j . The spatial string tension is chosen in the form 



./^jT)=c^g\T)T, (45) 

with the two- loop expression for g'^{T) 



-2/,N ^, 1 ^ , ^1 1 /^i M , _ T 



where 



^-^(t) = 26oln— + ^ln(^21n — j, t = - (46) 



bo={'iN.-lN,)' h=('4N!-{'4N.~ir)Nf 



3 3 V 167^2' V3 ' 3 N/ V (167r2)2- 

The measured in [26^ spatial string tension in pure glue QCD corresponds to the values 
of Ca = 0.566 ± 0.013 and L^ = K/T^ = 0.104 ± 0.009. On the left panel of Fig. El are 
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T/T, 



Figure 2: Left panel: The temperature over the square root of the spatial string tension 
versus T/T^ for pure glue QCD. Solid line corresponds to Eq.(p3|). The lattice data are 
from [26]. Right panel: Chromoelectric Debye mass mn/T for 2-fiavor QCD (upper lines) 
and quenched {Nf = 0) QCD (lower lines) versus T/Tc. Solid lines are calculated using 
muiT) = 2.09 ^yaJT), where ■\/crs(T) corresponds to Eq. pH|l with Nf = 2 for upper 
solid line and Nf = for lower solid line. Dashed lines are calculated using Eq. ljTfj) . 
Nf = 2-upper line, Nf = 0-lower line. The lattice data are from [13J. 



shown lattice data [2^1 and the theoretical curve (solid line) for T/ y'asiT) calculated 
according to ^ with c^ = 0.564, L^ = 0.104 and Nf = 0. 

On the right panel of Fig. |21 are shown lattice data [13j and theoretical curves for 
the Debye mass in quenched {Nf = 0) and 2-flavor QCD. Solid lines correspond to our 
theoretical prediction, nij:,{T) = co^/o^iT)-, with cc = 2.09 and for ^Jog (T) we exploit 
the same parameter (c^ = 0.564, L^- = 0.104) as in the left panel. The upper solid line 
is for the Debye mass in 2-flavor QCD, and the lower - for quenched QCD. We note 
that in computing mD(T) using (^3]) . P^l) all dependence on Nf enters only through the 
Gell-Mann-Low coefficients Bq and bi. For comparison we display in the right panel of 
Fig. 121 dashed lines for mDiT)/T, calculated with a perturbative inspired ansatz [15] 



m 



Latt 
D 



{T)=A 



N, 



Nf 
-f9iT)T. 



-Nf=0 



(47) 
1/(1.14 ■ 27r) Jig, and for the 



Quenched QCD corresponds to A^^^q = 1.51 and La^~ — 

2-flavor QCD Ajv^=2 = 1-42 and L^^^^ = 1/(0.77 ■ 27r) ^. 

Let us now consider higher orders of BPTh for ttid- From the gauge-invariant ex- 
pansion Q one obtains the next term {W^^\r, P)), which contains the double gluon 
propagator {afj,-^{xi)a^2{x2)a,^j^{yi)a^2{y2))a in the background field of the Wilson loop, 
which is proportional to g^{T). One can show that the background averaging of this 
propagator attached to the Wilson loop yields the diagram of the exchange of a double 



"'^ Physical justification for resorting to diniensionally reduced regime at T = ^/ajT) was given by 
(see also |2Z1) 
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gluon gluelump between Q and Q, and therefore the NLO BPTh Debye mass will coin- 
cide with the double gluon gluelump mass, computed for T = in [31 analytically and in 
PH on the lattice. As a result the lightest 2g gluelump mass appeared to be 1.75 times 
heavier than the lightest Ig gluelump mass. We expect therefore that also at T > Tc the 
same ratio of masses takes place, so that the asymptotics of gluon (gauge-invariant and 
background-averaged) gluon exchanges in BPTh has the form [^ EH! 

yGS(r, T) = A'^e-<''^^'^> - £aMi^i!l!e— i,'^"C^)- + . . . , (48) 

' 3 r r 

where C2(r) contains the asymptotic freedom logarithm. 

One can see in ()48p that the second term on the r.h.s. is subleading and small as 



compared to the first one, both due to {os Y and due to higher mass of m\f . At 
this point it is essential to note that this second term should enter as a sum over all 
possible 2g gluelumps. In one particular case, when two gluons form a color singlet, they 
decouple from the plane surface of the Wilson loop and create a 2g glueball, coupled by 
the spatial string (see FigOl). The corresponding glueball mass is computed in [22] and is 
1.7 times larger than the LO BPTh mass (j44|) . which is denoted as tti)-,^ (T) in (|48p . It 
is interesting that the aforementioned glueball mass corresponds to the gauge-invariant 
Debye mass suggested in [10], and as seen in (jIHj) it appears in the NLO BPTh, giving 
a small correction to the LO Debye screening potential. Therefore one can identify the 
Debye mass mr) = nijf with the accuracy 0{as 'e~^'^^'^) where Aitid = rn\)^ ~m\)^ > 
O.QGeV. 




Figure 3: NLO BPTh contribution to the Debye screening potential from the 2g glueball 
exchange 

The consideration above was done for the chromoelectric Debye mass, which appears 
in the screening Coulomb potential in the temporal plane (4i) appearing due to the 
(j44(x, y) gluelump Green's function. One can similarly consider exchange of "magnetic 
gluon" by insertion of the magnetic field vertex Fik ~ -Diflfc — Dj-ai into the Wilson or 
Polyakov loops. This vertex automatically appears in the Green's function from the 
term oc exp^ga^^i, J^ F^^dr) creating spin-dependent interaction. The same procedure as 
above leads to the "magnetic gluelump" Green's function, which differs from the "electric 
gluelump" case by the nonzero gluelump momentum L = 1. The corresponding mass is 
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easily obtained as in (P3|). giving 



£l(/^o, ^o) = 2y (1 + L + 2nr)aadj, £i = ^2^0 ^ 4:y/a's. 
Thus nonperturbative magnetic Debye mass is \/2 times heavier than the electric one. 

5 Conclusions 

We have studied Debye screening in the hot nonabelian theory. For that purpose the 
gauge-invariant definition of the free energy of the static QQ-pair in the singlet color state 
was given in terms of the thermal Wilson loop. Due to the chromomagnetic confinement 
persisting at all temperatures T, the hot QCD is essentially nonperturbative. To account 
for this fact in a gauge-invariant way the BPTh was developed for the thermal Wilson 
loop using path-integral FFS formalism. As a result one obtains from the thermal Wilson 
loop the screened Coulomb potential with the screening mass corresponding to the lowest 
gluelamp mass. Applying the Hamiltonian formalism to the BPTh Green's functions with 
the einbein technic the gluelump mass spectrum was obtained. As a result, we have derived 
the leading term of the BPTh for the Debye mass which is the purely nonperturbative, 
mD{T) = cd\/oJT) with cd ~ 2.06. 

Comparison of our theoretical prediction (solid lines on the right panel Fig. 12)) with the 
perturbative-like ansatz ()47|) (dashed lines) shows that both agree reasonably with lattice 
data in the temperature interval Tc < T < 5Tc] the agreement is slightly better for our 
results. At the same time, in (j77j) a fitting constant is used A^f ~ 1.5, which is necessary 
even at T/Tc ~ 5. At this point one can discuss the accuracy and approximations of 
our approach. As it was checked in numerous applications to hadron masses and wave 
function (see review ^7j) the accuracy of the Hamiltonian technic is around (5 -^ 7)%, 
while the area law is as accurate for loop sizes beyond A^ ~ 0.2 fm. At smaller distances 
the area law in ([THjl . (|T!Hl is replaced by the "area squared" expression [inj which yields 
effectively much smaller m£)(T). Therefore we expect that the Debye regime (^ with 
the rriD as in ()44j) starts at r > Ag ~ 0.2 fm. As a whole we expect the accuracy of the 
first approximation of our approach, Eq. ()44|) to be better than 10%, taking also into 
account the bias in the definition of a'f for the gluelump. The temperature region near 
Tc needs additional care because i) the behaviour (jl^ deviates from the data (see left 
panel of Fig. ^ and ii) contribution of chromoelectric fields above T^ (correlator Df, see 
|27j ) which was neglected above. Both points can be cured and will be given elsewhere. 
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